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1. Motivation 

It has been speculated that the observation of strange quark matter droplets might 
be an unambiguous way to follow up the transient existence of a quark gluon plasma 
(QGP) expected to be created in ultrarelativistic heavy ion collisions. Due to the Pauli 
principle, replacing several up and down quarks (u,d) by strange quarks (s), will lower 
the Fermi energy and the mass of strangelets consisting of multi-quark droplets. Their 
mass can thus become smaller than the mass of a corresponding ordinary nucleus (with 
the same baryon number B) and may then correspond to a stable state |^, |^. The 
mass per baryon of a strangelet may also be intermediate between the A mass and the 
nucleon mass thus producing a metastable state which then cannot decay into A’s 0. 
Both predictions are found in the MIT bag model, when its energy density constant 
is varied {B is the crucial parameter). Large values (> 180MeV) lead to 
large quark matter droplet masses and therefore to instability versus strong decays. 
Strangelets are found to be stable for B^^'^ < 150 MeV and metastable for intermediate 
values [g |. 

Several experiments are carried out at the Brookhaven AGS (E864, E878) and 
the GERN SPS (NA52), which look for strangelet production (small Z/A ratios) 
P, 0, |[. Here we present a novel method, which can allow for the measurement of 
the transient strange quark matter state, even if it decays on strong interaction time 
scales. The method is based on unequal particle correlations; they are sensitive to the 
time dependent charge and space-time expansion of the system. 


2. The Two-Phase Thermodynamical Model 


The dynamical evolution of the mixed phase consisting of a quark gluon plasma and 
hadronic gas can be described in a two-phase model which takes into account equilibrium 
as well as non-equilibrium features [|, p!0f| . 

Two main assumptions are made : 

1. The QGP is surrounded by a layer of hadron gas and equilibrium is maintained 
during the evolution. 

The Gibbs equilibrium conditions write: 
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2. Non-equilibrium evaporation is incorporated by a time dependent emission of hadrons 
from the surface of the hadronic hreball. The rate of frozen-out hadrons is assumed to 
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be governed by the actual hadronic time dependent densities with an universal effective 
relative rate F [p. The volume change (did) during the evolution of the plasma is 
determined by the thermodynamical energy according to: 


dE = TdS -pdV + 


( 2 ) 


In this model the volume increase of the expanding system competes with the volume 
decrease caused by freeze-out emission. Within this model, it is possible to follow the 
evolution of mass, entropy and strangeness fraction of the system. From this, relative 
yields of particles can be extracted which compare reasonably well with experimental 
data 1]^ . 

Fig.l shows the predicted evolution of the quark chemical potential /iq. It decreases 
from /iq 115 MeV to ~ 15MeV during the hadronization process. Conversely, 
the s quark chemical potential /is increases from /is = 0 MeV to several tens of MeV. 
Consequently, the strangeness fraction fg = (W — Ns)/Atot increases. This leads to 
a strange multi-quark object if the bag constant jg sufficiently small to allow for 
cooling of the system (see fig.2). 

An important new aspect of these calculations is the fact that the system necessarily 
leaves the usual /iq —T plane of the phase diagram. It enters quickly into the strangeness 
sector. Two important features emerge from these calculations (see fig.3 as an 
example): 

1) Strange and anti-strange quarks, produced in equal amount in the hot plasma, do 
not hadronize at the same time (in a baryon rich system). This so-called ’’distillation 
process” [O] results in the formation of stable or metastable blobs of strange matter in 


the case of low bag model constants. (In this case, a cooling of the system is predicted.) 
The evaporation from the surface of the hadron gas which is rich of anti-strangeness 
carries away etc. thus charging up the remaining mixed system with net positive 
strangeness. 

2) The large fg (~ 1.5 — 2) attained imply negative charge-to-mass ratios {Z/A = —0.1 
for stable, ZIA = —0.45 for metastable droplets). 

In the present work, these time separation and negative charge values characterize 
the transient existence of the plasma. They can be used to search for the distillation 
process by correlations of unlike particles. 

3. Influence of Coulomb Interactions on Particle Interferometry 

In the case of charged particles (like , K~), their Coulomb interaction modifies the 
plane and outgoing scattered waves and the two-particle amplitude takes the form: 


= e*'y/IM fe-*^^”F(-*r/, 1, tp) + $^(r) 


( 3 ) 
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where 6 is the Coulomb s-wave phase shift, F{—iri, 1, i^) is the confluent hypergeometric 
function, p = kr + kr, rj = {ka)~^, |a| is the Bohr radius of the two-particle system, 
<h^(r) is the scattered wave, and 

Aij]) = 27r?7[exp(27rr7) — 1]“^ (4) 


is the well known Coulomb factor. The A{ri) factor differs from unity mainly for 
\2Tir]\~^ < 1 which means fc < 12 MeV/c in the case of the K^K~ system (a = —110 
fm). At high energies the two-body Coulomb effect is often treated through a size- 
independent correction factor A{ri) to the usual quantum statistical correlations (i.e., it 
is assumed that r /\a\ -C 1 so that F{—iri, I, ip) ~ 1). It is thus not used to study the 
space-time evolution of the system. Contrary, at low energies this effect becomes the 
main tool for the study of this evolution due to large transit times. In addition, the 
charged source itself can also affect the above two-particle correlation pattern due to the 
long range interaction (e.g. accelerating and deaccelerating K~ for Z > 0). This 
effect is much more clearly visible in the single particle spectra, while it nearly cancels 
out (by taking the denominator) in the correlation function. 

Experimental evidence and detailed analysis of such Coulomb two- and three-body 
effects on two unlike particle correlation functions can be found in |T3 . The three-body 


effect itself was observed experimentally more than ten years ago [HI in intermediate 
energy heavy ion collisions and quantitatively analyzed in reactions producing light 
particles or unstable fragments (Li*, B*) and short-lived resonant states. If their 
decay occurs into two different particles (with different charge/mass ratios) within the 
Coulomb held of the residual nuclear system (target/projectile fragments), this held can 
substantially distort their correlation function. Thus a complete three-body calculation 
is required to account for the ehect of hnal state interactions. One way to look at this 
width/lifetime ehect is to compare correlation functions measured for diherent relative 
velocities Bin- One then observes a shift of the resonance peak which rehects the 
diherent Z/A values of the coincident particles which suher the Coulomb accelerations. 
This ehect has been later used to explain a variety of two particle correlation patterns 
in the intermediate energy domain [TO. In the adiabatic limit, it has also been studied 


within a detailed and full three-body quantum approach HH- This ehect is found 


to decrease for particles emitted at large relative space-time intervals or with large 
momenta with respect to the source. It is increasing with increasing source charge and 
also with increasing charges and masses of the particles. 


4. Access to Time Sequence of Emission nsing Unlike-Particle 
Correlations 

In addition to the known directional dependence of the correlations of two identical 
particles, it has been shown (at low energy heavy ion collisions) that the hnal state 
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interactions between non-identical particles can provide information about the time- 
order of emission. Thus, based on classical trajectory calculations, it was proposed to use 
the velocity difference spectrum between two different charged particles and examine its 
dependence on particle energy in conjunction with the corresponding relative momentum 
correlation function |]^. Independently, in the quantum approach [^, ^ it was 
shown that the anisotropy of the space-time distribution is reflected in the directional 
dependence of unlike-particle correlations and can be directly used to measure the delays 
in the emission of particles of different types. 

With the assumption of a sufficiently small density in momentum space, the 
correlation is only due to the mutual hnal state interaction; the two-particle amplitude. 


in the absence of the Coulomb interaction, takes the form pi 


vy = e 


—ikr 


+ *hfc(r), 


( 5 ) 


where k = q/2 and f are the relative momentum and spatial coordinates of the two 
particles in their cms system, <hfc(r) is the scattered wave. The correlation function for 
two non-identical particles is 

1 + R{Pi,P 2 ) = 1-1- (|<hfc(r)p-|-2 Re<hfc(r) cos(fcr) — 2 Im <hfc(r) sin(fcr)) .(6) 


The directional dependence in the correlation function is contained in the odd term, 
Imd>fc(r) sin(fcrj). In the limit of large relative emission times t = ti — t 2 , u|t| 3> r, the 
vector r is only slightly modified by averaging over the spatial location of the emission 
points in the rest frame of the source. So the vector f is nearly parallel or antiparallel 
to the pair velocity v depending on the sign of the time difference. The odd part of the 
correlation function is sensitive to this sign. This mean relative emission time, including 
its sign can be determined by comparing the correlation functions R'^ and R~ [^|. In 


these functions the sign of the scalar product k ■ v is hxed to be larger or smaller than 
zero (using their relative angle) respectively. Note that for particles of equal masses the 
sign of A; -u is almost equal to the sign of the velocity difference Vi—V 2 - This means that 
the interaction of the two particles will be different in the case where the faster particle 
is emitted earlier as compared to the case of its later emission. 

In the case of unlike charged particles, additional odd terms in the k ■ r product 
appear due to the confluent hypergeometrical function. This modihes the plane wave 
and preserves the sensitivity to the time difference even in the case of a weak effect of 


the strong interaction. It has also been shown |^| that the Coulomb field of the residual 
nucleus slightly affects this result. 

These predictions have been compared to experimental (p,d) correlations measured 
for the Nb-|-Pb, Xe-fTi, Xe-fSn collisions studied at GANIL [^, |^, The 

ratio (1 -|- R'^)/{1 + R~) deviates from unity at low q. This can be attributed to 
an earlier emission of deuterons rather than protons on the average. Thus, it has 
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been experimentally demonstrated for the first time that the study of the directional 
asymmetries in the correlations of non-identical particles provides a unique tool to 
measure the time-ordering of particle emission. This technique will here be extended to 
emission time differences between strange, anti-strange and non-strange particles . 
We study its sensitivity to the time delays related to the production of the transient 
strange quark matter state. 


5. Calculations with an Event Generator 


A quantum two-particle final state interaction code is used together with a three- 
body calculation m to simulate the possible influence on K^K correlations caused 
by 


1) different distributions of the and K~ emission times and 

2) by different charge states of the emitting source. 

Results are presented for an event generator describing a mid-rapidity source with 
transverse mass spectra conditions typical for AGS experiments. The source is assumed 
to be Gaussian in space coordinates and exponential in time. Different sets of time 
origins are used for the two particles. 

Fig.4 shows some simulations with mean emission time differences between and 
K~ of —5 to -l-lOfm/c. The figure shows the ratio of the relative momentum correlation 
functions (1-|-R ^)/{1 + R~) as obtained with the above selections for values of the scalar 
product fc ■ r > 0 and < 0, respectively. Glearly, ratios less than one are associated with 
a scenario where the K~'’s are emitted later on the average. 

Fig.5 shows the dependence of the above ratio on the charge of the source. The 
time difference is fixed to -l-lOfm/c. The upper diagram shows the result of a 2-body 
calculation. The two lower diagrams correspond to 3-body calculations with source 
charges Z = 50 and Z = —50, respectively. We have checked that the neglection of 
the strong hnal state interaction only slightly modifies these results. Therefore, most 
of the effect originates from the Goulomb interaction between the two kaons with an 
additional slight enhancement (suppression) of this deviation from unity created by a 
positive (negative) source charge. 


6. Calculations with the Two Phase Thermodynamical Model coupled 
to Final State Interactions 

The previous section clearly shows the domain of sensitivity of the K^K~ correlations 
to time differences and charge effects. In order to do similar calculations with a 
more sophisticated particle source, we have described the source evolution using the 
thermodynamical model presented in section 2. Three sets of parameters have been 
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chosen: 

set 1 h 2 (representing Au+Au at AGS energies): 

initial mass A = 394, entropy per baryon S/A = 10, initial net strangeness fs = 0; 

set 1: = 160 MeV, strangelets (albeit metastable) are formed. 

set 2: = 235 MeV, strangelets are not distilled, their thermal formation is also 

highly suppressed. They are not (meta-)stable. 

set 3 ( representing S+Au at SPS energies): 

bag model constant = 235 MeV, 

initial mass A = 100, entropy per baryon S'/A = 45 , initial net strangeness fs = 0. 
Corresponding calculated time evolutions are shown in hgs. 6 and 7 for particle yields, 
initial quark blob masses and charges. In the case of a high bag constant the system 
heats up slightly. The entropy per baryon in the QGP is larger than in the hadronic 
phase. This results in a fast hadronization (see the radius and charge time evolutions). 

are emitted earlier than K~. For low bag constants, a cold strangelet emerges in 
a few tens of fm/c with a radius of ~ 2.5 fm, a hnite net strangeness content and a 
negative charge. 

These predicted time dependent yields are used as input for the hnal state interaction 
calculations. Fig.8 shows the calculated correlation function ratios (1 + R^)/{1 + R~) 
for K^K~ pairs. Similar functional dependencies are found as with the simple event 
generator calculations. 

The ratio of the measured correlation functions for k ■ v > t) and k ■ v < t) shows 
that the K~ is emitted later in the case of a low bag constant value. The effect is seen 
as a 30 — 35% dip relative to equal mean time production (i? = 1). It extends up to 
q K, 30MeV/c for the parameters of set 1. 

A weaker effect (10%) is found, extending up to g ~ 40 MeV/c when set 2 (high bag 
constant) is used. The effect is nearly gone for the higher initial entropy set 3. This may 
be attributed to freeze-out time distributions of and K~ which are quite similar as 
seen in the predicted yields (hgs. 3 and 7). The slight difference between the sets 2&3 
has to be attributed to the smaller Goulomb inhuence of the source itself (A = 100) 
in the entropy S'/A = 45 calculation. Our simulations (see hg.5 and discussion above ) 
support this explanation (including the direction of this effect) as long as an averaged 
positive charge is present (in agreement with the present model calculation displayed in 
fig-6). 

7. CONCLUSION AND PERSPECTIVES 

We have presented a novel method, which can allow for the measurement of the transient 
strange quark matter state, even if it decays on strong interaction time scales. The 
method is based on unequal particle correlations. They are sensitive to the time 



dependent charge and space-time expansion of the system. The dynamical evolution 
model for the mixed QGP-hadronic phase has been applied to describe the strangeness 
distillation process. This predicts a time dependent freeze-out and the occurence of low 
charge-to-mass ratios. The existence of strangelets can thus be characterized by the 
correlation approach of unlike particles presented here. 

Hydrodynamical calculations including the expansion of the system are currently 
underway. Also microscopic analysis of this effect, based on the UrQMD-simulations, 
are presently under study. If combined with the identical particle interferometry (HBT), 
this technique offers the possibility to study different expansion processes associated with 
different transverse motion and sizes. Correlations of other particles (pA, pvr”, etc.) are 
presently studied. 
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Figure captions 


Figure 1. Time evolution of the quark and strange chemical potential for parameter 
set 3. Hadronic and quark matter rapidly moving out of the (/iq,r) plane, entering 
the strange sector (fg ^ 0). The fit of a static thermal source which yields the same 
hadron ratios as the dynamical model is also depicted. 


Figure 2. Time evolution of baryon number strangeness fraction jg and 
temperature T for S'/Hinit = 25 and fg = 0.25 for a (a) high = 235 MeV) 

and (b) low = 145 MeV) bag constant, respectively. 


Figure 3. Yields for 7r+, K'^, K , p and A particles as a function of time for the AGS 
parameter sets 1&2 (upper figure B^/'^ = 160MeV, lower figure B^/'^ = 235 MeV). 


Figure 4. Simulation of different time origins for K'^ and K emission with the event 
generator. Shown are the ratios (1 + R^)/{1 + R~) as explained in the text. 
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Figure 5. Simulation of the influence of the source charge state on the ratio 
(1 + i?+)/(l + R~) with the event generator (see text). 


Figure 6. Radius and charge of the remaining quark matter blob as a function of 
time t as calculated with the two-phase thermodynamical model. 


Figure 7. Particle yields as a function of time t for a bag constant = 235 MeV, 
initial entropy S/A = 45, initial baryon number ^init = 100, representing an SPS 
parameter set (set 3 in the text). 


Figure 8. Ratios of correlation functions (Id- + R ) for three different 

parameter sets (see text). 
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